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Liquid Xenon (LXe) is an excellent material for experiments designed to detect dark matter in the 
form of Weakly Interacting Massive Particles (WIMPs). A low energy detection threshold is essential 
for a sensitive WIMP search. The understanding of the relative scintillation efficiency (Xeff) and 
ionization yield of low energy nuclear recoils in LXe is limited for energies below 10 keV. In this paper, 
we present new measurements that extend the energy down to 4 keV, finding that Xeff decreases with 
decreasing energy. We also measure the quenching of scintillation efficiency due to the electric field 
in LXe, finding no significant field dependence. 

PACS numbers: 95.35.+d, 29.40.Mc, 95.55.Vj 



I. INTRODUCTION 

Liquid xenon is increasingly used as the detection ma- 
terial in direct searches for WIMP dark matter [1] . Recent 
developments in tw^o-phase (gas/liquid) xenon detec- 
tors [2 -4 J has resulted in stringent limits on the WIMP- 
nucleon cross-section, constraining theories of physics 
beyond the standard model, such as supersymmetry. 
WIMPs will deposit a small amount of energy in the 
LXe through elastic scatters with xenon nuclei. Part of 
the deposited energy is converted into observable sig- 
nals of scintillation light and ionization electrons. The 
rest of the energy is converted into heat and can not be 
easily measured. Understanding these effects will help 
determine nuclear recoil energies and ultimately play a 
part in determining the WIMP-nucleon cross-section. 

In a two-phase xenon detector, two signals are mea- 
sured. The first is the direct scintillation light, denoted 
as SI. The second is the proportional scintillation light 
in the gas phase, denoted as S2, which is proportional to 
the ionization electrons that survive electron-ion recom- 
bination and are extracted into the gas. Figure [l] gives 
an illustration of the signal production and collection in 
a two-phase xenon detector. 

For a given event in the LXe, the nuclear recoil energy 
can be determined based on the scintillation signal SI 
ISO. However, it is much more convenient to calibrate 
the detector using electron recoil events. The tradition in 
the field |5 - 9 1 is to base the energy calibration on 122 keV 
electron recoils from a ^""Co gamma source. The relative 
scintillation efficiency, Xeff/ defined as the ratio between the 
electron equivalent energy (E^e) and the true nuclear re- 



coil energy becomes necessary for determining the 
nuclear energy scale and, therefore, the WIMP detection 
sensitivity. E^e is inferred from the scintillation signal 
yield due to monoenergetic electron recoils. Xeff has no 
units and is defined at zero electric field in LXe relative 
to 122 keV gamma rays. 

If an electric field is applied to the LXe, the scintil- 
lation yields for both electron and nuclear recoils are 
suppressed by additional factors S^ and S„, respectively. 
The relative scintillation efficiency can be calculated as 



Xeff — E^^/Ef • Sg/Sfi 



(1) 



The quantity Se for 122 keV electron recoils from a ^""Co 
source has been measured very accurately |10|. S„ has 
been measured for 56 keV nuclear recoils, with electric 
fields up to a few kV/cm in LXe [7, 10 1, but no measure- 
ment is available for nuclear recoils at other energies. 
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FIG. 1: (Color online) Illustration of the signal production and 
collection in a two-phase xenon detector. 
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Two methods have been utilized to determine £eff as 
a function of energy: a) Using a fixed-energy neutron 
beam experiment, detecting neutrons that scatter in the 
LXe at a known scattering angle, and b) Comparing neu- 
tron calibration data to Monte Carlo simulations without 
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tagging the scattered neutron. Using method a), Xeff has 
been measured by a number of groups for nuclear recoils 
above 10 keV ||5}Q. There are also two measurements of 
this type reporting results below 10 keV, one suggesting 
an increasing Xeff with decreasing energies |8|, and an- 
other indicating a roughly constant Xeff of 0.15 down to 
5 keV |9|. The XENONIO and the ZEPLIN-III collabora- 
tions have also determined Xeff with method b) [4 . 11 . 12J . 
In the XENONIO analysis, Xeff does not decrease much 
at low energies, while in the ZEPLIN-III measurement, 
the data imply a precipitous drop at low energies. 

In this paper, we report on measurements of Xeff at 
zero field and Sn at two different fields (0.73 kV/cm and 
1.5 kV/cm) for nuclear recoils between 4 and 66 keV in 
a single phase detector (Sl-only). We repeated these 
measurements using a dual phase detector (SI and S2 
signals) at 1 kV/cm in the liquid and 10 kV/cm in the 
gas as well as 4 kV/cm in the liquid and 8 kV/cm in the 
gas. With the dual phase detector we also measured the 
ionization signal yield for nuclear recoils. We present the 
experimental apparatus in Sectionlllj the data analysis in 
Section [In) the results in Section [iVf a theoretical model 
of Xeff in Section |V| and a discussion of the results in 
SectionlVll 



II. EXPERIMENTAL APPARATUS 

The measurement was performed with a setup com- 
prising a deuterium-deuterium neutron generator 1 13 1, a 
LXe detector, and an organic liquid scintillator detector, 
as shown in Figure |2] The neutron generator produces 
2.8 MeV neutrons at a rate of 10^ n/s. The liquid scin- 
tillator detector is a BC501A organic scintillator module 
3.8 cm in diameter and 3.8 cm in height, viewed by a pho- 
tomultiplier (PMT). Both the neutron generator and the 
organic scintillator detector have previously been used to 
perform measurements of nuclear recoils in liquid argon 
[14J and liquid neon [15]. The LXe detector is made of a 
cylinder of LXe viewed by two Hamamatsu R9869 PMTs, 
as shown in Figure |3] The PMTs are specially designed 
for LXe applications. They have a bialkali photo-cathode 
and a quartz window with an aluminum strip pattern on 
the photo-cathode. The two PMTs have a quantum effi- 
ciency of 36% for LXe scintillation light at 175 nm. The 
collection efficiency from the photo-cathode to the first 
dynode is about 70%. The active LXe target is 5 cm in 
diameter and 2 cm in height and is surrounded by poly- 
tetrafluoroethylene (PTFE) for UV light reflection. The 
thickness of the PTFE is minimized (11.5 mm) to reduce 
neutron multiple scatters with surrounding materials. 
Two stainless steel mesh grids, each with 90% optical 
transparency, are installed to apply electric field to the 
LXe. In dual phase mode, a third grid is added to apply 
a separate electric field in the xenon gas region. 

The LXe detector is located in an aluminum vac- 
uum cryostat and cooled by a pulse-tube refrigera- 
tor (PTR). The cryogenic system is described else- 



where fT6l. The neutron generator is shielded by 
30.5 cm X 30.5 cm x 30.5 cm water boxes and two 10 cm 
X 30 cm X 30 cm polyethylene slabs to block and ab- 
sorb neutrons that are emitted in directions other than 
toward the cryostat. The distance from the neutron gen- 
erator to the center of the LXe was 76 cm. The distance 
between the center of the LXe and the center of the or- 
ganic scintillator varied from 16 to 20 cm between runs. 
The scattering angle, defined by the position of the or- 
ganic scintillator (Figure |2]), was varied from 25 to 125 
degrees to vary the associated recoil energy. 
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FIG. 2: (Color online) The setup for nuclear recoil scintillation 
efficiency measurement in LXe. Not drawn to scale. 
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FIG. 3: (Color online) Schematic of the dual phase LXe detector. 
Spaces not drawn were filled with PTFE pieces. PMTs, LXe and 
xenon gas regions drawn to scale. 

During the single phase runs, the PMT waveforms 
were recorded with an 8-bit oscilloscope, model TDS- 
5034B from Tektronix. The oscilloscope's logic gate was 
used to trigger the data acquisition system, triggering 
on triple coincidence (both SI signals and the organic 
scintillator signal) in a 80 ns window. 
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For the dual phase runs, a VME 12-bit, 250 MS/s dig- 
itizer (CAEN VI 720) was used, because of its higher 
dynamic range. An external trigger for the VME digi- 
tizer was generated using external NIM modules. In this 
mode the data acquisition system was triggered by the 
S2 signals. At a few scattering angles, data acquisition 
was instead triggered by the SI signals to test for any 
systematic effects due to the trigger. For the S2 trigger, 
shown in Figurelij the PMT signals were added and inte- 
grated using a FAN IN/OUT and an integrator module. 
The summed and integrated signal then went to a dis- 
criminator to select the S2 signals. The scintillator signal 
was sent through a discriminator to avoid small pulses. 
The output signal was used to generate a 15 /is pulse with 
the help of a gate generator. This pulse in coincidence 
with the S2 pulse generated the trigger for the acquisi- 
tion system. For the SI trigger we reproduced the trigger 
system described above for single phase operation, but 
using NIM modules. In both the single and dual phase 
runs, the data were recorded at a rate of ~5 Hz. 



the light yield drops to 9.5 + 0.2 pe/keVee at zero field 
and 4.3 + 0.1 pe/keVee at 1.0 kV/cm drift field. The light 
yields from the 122 keV events were monitored over the 
entire period of the measurements finding an additional 
fluctuation less than 3%. 
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FIG. 4: (Color online) Trigger system used in the dual phase 
runs. The PMT signals are summed and integrated to select 
the S2 signals. The organic scintillator signal generates a 15 
/is pulse. This pulse and the S2 signals in coincidence trigger 
the acquisition system. For the single phase runs, triple co- 
incidence of the SI signals and the organic scintillator signal 
triggers the acquisition system. 



FIG. 5: (Color online) Scintillation light yield for 122 keV 
gamma rays in the LXe detector running in single phase. The 
left histogram was taken at 1.5 kV/cm yielding 4.8+0.1 pe/keVee 
and a resolution (o/E) of 18%. The right histogram was taken 
at 0.0 kV/cm yielding 10.8 + 0.1 pe/keVee and a resolution of 
8.8%. 



During the single phase runs, the LXe purity was mon- 
itored by measuring the stability of the scintillation sig- 
nal yield from 122 keV gamma rays. During the dual 
phase runs, the purity was also monitored by measuring 
the electron lifetime, t, found by fitting the S2 spectrum 
with S2 = S2o exp [-dt/z] where dt is the electron drift 
time, measured as the time between the SI and S2 sig- 
nals. For the data presented here, electron lifetime was 
greater than 40 jis and continuously improved over the 
course of the experiment. By the end of the experiment, 
the electron lifetime was 90 /is. The total drift time for 
events at the cathode grid is 12 /is. The S2 signals were 
corrected for the electron lifetime. Figure [6] shows a typ- 
ical S2 spectrum for 122 keV gamma rays taken during 
the dual phase runs. 



Throughout the runs, periodic calibrations were per- 
formed to test the stability of the PMTs and measure the 
purity of the LXe. The gains of the two PMTs were 
measured from the single photoelectron (pe) spectra 
by using light emitted from a blue LED located inside 
the LXe detector. The energy scale is calibrated using 
122 keV gamma rays from a ^^Co source located out- 
side the cryostat. The scintillation signal yields for the 
122 keV gamma rays in LXe in the single phase detec- 
tor are 10.8 + 0.1 pe/keVee (keV in electron equivalent) 
at zero field and 4.8 + 0.1 pe/keVee at 1.5 kV/cm (Fig- 
ure |5|. After adding a third grid and PTFE spacer and 
then removing some LXe to run the dual phase mode. 



in. FINDING Xeff 

To obtain the Xeff value for each experimental setup we 
take the following steps. First, apply a set of cuts to re- 
duce uninteresting events such as noise events, gamma 
ray scatters, neutron inelastic scatters and multiple elas- 
tic scatters. Next, the energy spectrum is obtained based 
on the light yield measurements and fitted to a spec- 
trum predicted through Monte Carlo simulation. The 
following subsections explain these steps. 
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FIG. 6: S2 spectrum for 122 keV gamma rays in the dual phase 
detector. 



A. Data Analysis 

To remove uninteresting events two sets of cuts were 
used. The first set removes noise events and events 
outside the energy window of interest. The second set 
consists of two cuts used to select single elastic nuclear 
recoil events (Figure [7]): 

• The first cut is based on the pulse shape discrim- 
ination (PSD) of the organic scintillator. The PSD 
is based on the relation between pulse height and 
pulse area. This is an effective way to separate the 
neutron events from gamma events in the organic 
scintillator. 



• The second cut uses the fact that neutrons take a 
longer time to travel from the LXe detector to the 
organic scintillator than gamma rays. The cut uses 
the time of flight (ToF) to remove events triggered 
by gamma rays or accidental coincidence. The cut 
selects the first half of the ToF peak because the 
contribution from scatters other than single elastic 
ones is negligible, as determined from Monte Carlo 
simulations. 



The nuclear recoil energy distributions (SI signals) af- 
ter these two cuts are shown in Figures [8] (a) and (b) for 6 
keVr nuclear recoils at zero electric fieldand 66 keVr nu- 
clear recoils at 1.5 kV/cm, respectively. Similarly, Figures 
|8](c) and (d) show S2 distributions for 6 keVr nuclear re- 
coils at zero electric field and 66 keVr nuclear recoils at 
1.5 kV/cm, after the cuts. In the 56 and 66 keVr runs it 
is easy to separate the background tail from the single 
elastic signal, as shown in Figure[8](b). This tail has been 
removed (SI < 50 pe) to generate the S2 spectrum in 
Figure [8] (d), and to find the Xeff values. 
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FIG. 7: (Color online) Two cuts used in the analysis to select 
single scatter nuclear recoil events, (a) pulse shape discrim- 
ination in the organic scintillator, (b) neutron time of flight 
between the LXe and the organic scintillator. Both plots show 
experimental data for 56 keVr run. 



B. Monte Carlo simulation 

The detector's response to neutrons is modeled using 
a Geant4 [17] simulation that takes into account the real- 
istic setup of the experimental apparatus as described in 
Sectionlin including the water shield around the neutron 
generator, the polyethylene shield around the organic 
scintillator detector, the aluminum cryostat, the stain- 
less steel cell, PTFE structure, PMTs, grids and LXe in 
the detector. The simulation stores the neutron scatter- 
ing position, time, energy and type of events that deposit 
energy in both the LXe detector and the liquid scintilla- 
tor. 

For the simulations we follow the work by fT2\ and 
used the Xe(n,n)Xe scattering cross-sections from the up- 
dated ENDF/B-VII data, instead of the ENDF/B-VI cross- 
section used by default with the software. As discussed 
in [18], a +3% uncertainty associated to the well depth 
parameter in the Optical Model Potential used to cal- 
culate the elastic cross-sections, leads to a conservative 
uncertainty which translates into an uncertainty up to 
+3% in the nuclear recoil spectra. The latter is included 
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FIG. 8: SI spectra in pe for (a): 6 keVj- nuclear recoils at zero electric field and (b): 66 keVj- nuclear recoils at 1.5 kV/cm, after 
applying the cuts. S2 spectra in pe for (c): 6 keVj- nuclear recoils at zero electric field and (d): 66 keVj- nuclear recoils at 1.5 kV/cm. 



when calculating the Xeff systematic uncertainty. 

Neutrons can deposit energy in LXe via elastic scatter- 
ing, inelastic scattering, or a mixture of both. For most of 
the events that satisfy the ToF cut between the LXe and 
the liquid scintillator, neutrons come directly from the 
neutron generator, make a pure single elastic scatter in 
LXe and reach the liquid scintillator. Single elastic events 
give a peak (see Figure |9l at the energy Er determined 
by the kinematics according to: 



Er ^ En- -^(1 - COS 6) 



(2) 



where En is the incoming neutron energy (2.8 MeV), 
while rUn and mxe are the masses of the neutron and 
the Xe nucleus, and 6 is the scattering angle. The spread 
of the energy deposition peak is from the spread of 6 due 
to the geometric width of the LXe detector, the width of 
the organic scintillator and the distance between them. 

Some of the neutrons from the generator may scatter 
first in other materials (e.g. PTFE) before entering the 
LXe detector. Neutrons may also scatter more than once 
in the LXe detector. These "multiple-scattering'' events 
have a variety of scattering angles and, therefore, de- 
posit a wider range of energies than the single elastic 



scatters. This contributes to the background tail under 
the pure single scattering peak. Through minimization 
of non-active material, the geometry of the detector was 
designed to reduce background from multiple scattering 
events, especially for runs performed at low scattering 
angles. Inelastic scattering events are very few and make 
negligible contribution to the energy spectrum (less than 
1%). A detailed description of the contributions of each 
background for each energy tested can be found in |19|. 

Although the LXe detector gives a high scintillation 
yield, detecting and resolving the peak for low energy 
nuclear recoils is still challenging. The energy spread for 
nuclear recoils below 10 keVr is dominated by statistical 
fluctuation of the photoelectrons in the PMTs. The trig- 
ger threshold, the trigger efficiency and the SI -finding 
algorithm efficiency must be taken into account in the 
analysis since not all events are effectively detected. In 
our measurement, the trigger of the LXe signal is from 
coincidence of the two PMTs with leading-edge discrim- 
inators. Effects of the trigger and software efficiencies 
are included in the Monte Carlo spectrum. A realis- 
tic model of the trigger efficiency was determined by a 
Monte Carlo simulation that included photon distribu- 
tion into the two PMTs, quantum efficiencies, statistical 
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FIG. 9: (Color online) Simulated distribution of energy deposi- 
tion in the LXe for 2.8 MeV neutrons that scatter and are tagged 
by the organic scintillator, for three different angles, assuming 
the same ToF cuts that were used in the experimental data. The 
solid line corresponds to all types of events, the blue dashed 
line to single elastic nuclear recoils, and the purple dotted line 
to multiple elastic scatters. Multiple scatter events and inelas- 
tic events are negligible compared to the single scatter events. 
The energy and uncertainty values indicated in the figure leg- 
ends are based on the mean and variance (one sigma) of the 
peak for all events. 



sampling of typical noise and single photoelectron w^ave- 
forms in the PMTs, the electronic trigger thresholds, and 
the SI -finding algorithm in the analysis software. Fig- 
urellOlshows the overall trigger efficiency when combin- 
ing tKese effects. 
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FIG. 10: (Color online) Detector triggering and software effi- 
ciencies for the single ( ) and dual phase (•) runs. The ■ 
points show the ideal case when the PMTs and the software 
have 100% efficiency for a single photoelectron. 

For each energy studied, the Xeff value is found by 
comparing the Monte Carlo generated spectrum with the 
measured spectrum, using a test according to equa- 
tion IM: 



X 



-L 



{ni-Vif 



(3) 



Z=l 



where N is the total number of bins, while Ui and Vi 
are the measured and Monte Carlo generated number 
of events in each energy bin. To perform the test, 
the total number of events in the Monte Carlo spectrum 
was normalized to that of the measured spectrum. The 
energy distribution (or Vi at different energy bins) from 
the Monte Carlo spectrum varies with different input of 
Xeff values. The best-fit Xeff value is obtained by min- 
imizing the x^ parameter. Figures 11 (a) and (b) show 
the data (points with error bars) and the Monte Carlo 
spectrum (line) after minimization, for 6 keVr nuclear 
recoils at zero electric field and 66 keVr nuclear recoils 
at 1.5 kV/cm. These plots correspond to the SI spectra 
shown in Figures [8] (a) and (b). Figures 11 (c) and (d) 
show the x^ vs. Xeff histograms. 



IV. RESULTS 

The largest uncertainties in Xeff derive from the x^ 
analysis and from the energy resolution applied to the 
Monte Carlo spectrum. The energy resolution is a com- 
bination of Poisson fluctuations in the PMTs, the PMT 
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FIG. 11: (Color online) SI measured spectrum (data points) and Monte Carlo comparisons (solid line) for (a) 6 keVj- nuclear 
recoils at zero electric field and (b) 66 keVj- nuclear recoils at 1.5 kV/cm. The background tail (SI < 50 pe) for the 66 keVj- run has 
been removed for the test. Also shown are the ^ vs Xeff histograms for (c) 6 keVj- nuclear recoil at zero electric field and (d) 
66 keVr nuclear recoils at 1.5 kV/cm. The points show the different Xeff values tested while the curve is a fit to the points used to 
find the Icr errors on the Xeff value. 



gain fluctuations, the geometry of the cell {ogeo), the op- 
tical properties of the materials and the intrinsic resolu- 
tion of the LXe. The total energy resolution, o, defined 
as the root mean squared of the terms mentioned, was 
measured to be cr = (3.2 + 0.4) V^? with N the number 
of photoelectrons for 56 and 66 keVr. This relationship 
was assumed for determining all L^a values. Because 
the geometry resolution is already included in the sim- 
ulation, the Monte Carlo data were convolved using an 

energy resolution o' = ^Ja^ - a^^g. The gain fluctuations 

were measured from the calibration runs, while the er- 
ror due to the optical properties was determined from 
an independent light simulation. After subtracting the 
geometry resolution, the overall uncertainty on a' is es- 
timated to be +1.0 ^fN, and this is propagated through 
to compute the systematic uncertainty in Xeff- The to- 
tal errors are shown in Figure [12] and given in Table [l] 
Figures and 13 compare our Xeff results with previ- 
ous analyses. Since the data show no significant nuclear 
quenching due to the electric field, £eff is computed 



using data from all of the different runs, regardless of 
electric field. 

An electric field applied to LXe will suppress electron- 
ion recombination and thus reduce the scintillation yield. 
This field induced quenching is significant for electronic 
recoils. For example, at 1 kV/cm, scintillation light yield 
from 122 keV gamma rays in LXe is reduced to 53% of its 
value at zero field [7J. For nuclear recoils, a near-unity 
value of Sn has been measured at 56 keVr [7,101. Here 
we measured the field induced quenching for nuclear 
recoils as low in energy as 4 keVr. No significant field 
quenching was observed for any energy or electric field. 
The average field induced quenching factor, S„, for 56 
keVr nuclear recoils at an electric field of 0.73 kV/cm is 
about 95% as given in Table [l| 

From the dual phase data we can determine the ioniza- 
tion yield (number of electrons escaping recombination 
per unit recoil energy). This number is determined from 
the S2 peak position for the nuclear recoils (Figure Isld) 
for example) and the number of photoelectrons per elec- 
tron determined from the ^""Co calibration runs. Figure 
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FIG. 12: (Color online) Scintillation efficiency for nuclear re- 
coils relative to that of 122 keV gamma rays in LXe at zero 
field, comparing this work (•) to previous measurements from 
Arneodo (a) |5|, Akimov (□) |6|, Aprile (ik) |7|, Chepel (O) |8| 
and Aprile (0)l 9 1. Also shown is the theoretical model (dashed 
line) explained in Section [V] which includes the Lindhard fac- 
tor, an electronic quenching due to bi-excitonic collisions and 
the effect of escaping electrons. 
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FIG. 14: (Color online) Ionization yield as a function of re- 
coil energy. Shown are the measured values in this work at 
1.00 kV/cm (A) and 4.00 kV/cm (•), along with previously 
measured values at 0.10 kV/cm (O), 0.27 kV/cm (□), 2.00 kV/cm 
(a) and 2.30 kV/cm(0) from 1 10 1, error bars omitted for clarity. 
Also shown are the ionization yields calculated by comparing 
the XENONIO nuclear recoil data and the Monte Carlo simu- 
lations 1 12 1 for single elastic recoils at 0.73 kV/cm, using two 
different methods (* and ik). 
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FIG. 13: (Color online) Scintillation efficiency for nuclear recoils 
measured in this work (•) and the theoretical model (dashed 
line) compared to the scintillation efficiency found from the 
neutron calibration data by the XENONIO (top shaded area) 
[12J and the ZEPLIN-III (bottom shaded area) 14J collabora- 
tions. 



14 shows the energy dependence of the ionization yield 
measured in this work for 1.0 kV/cm and 4.0 kV/cm, as 
well as previous measurements [10 J and the calculated 
values when comparing the XENONIO nuclear recoil 
data and Monte Carlo simulations (12]. The ionization 



yield errors shown in Figure 14 were derived from the 
width of the S2 signals from nuclear recoils and from 
the ^""Co calibrations. By comparing the dual phase runs 
triggered by the SI signals with the runs triggered by 
the S2 signals, we determined that there is no significant 
uncertainty in the S2 signals due to the trigger. 



V. EMPIRICAL MODEL OF Xeff 

The data shown above reveal a relative scintillation 
efficiency that decreases with decreasing energy. A suit- 
able theoretical expression for Xeff in LXe can be written 
as the product of at least three components: 



Xeff — ^ncl ' ^esc ' ^el 



(4) 



First is the Lindhard factor [21], q^xch which quantifies 
the larger fraction of energy dissipated into atomic mo- 
tion or heat in a nuclear recoil compared to that for an 
electron recoil. As a function of recoil energy, E^/ the 
Lindhard factor can be written as 



^?ncl 



(5) 



where for a nucleus with atomic number Z and mass 
number A, fc = g{e) = S.Oe^-^^+O.Ve^-^+e, 

where e is the reduced energy e = 11.5 • £r • Z~^^^. 

The second term, ijesc/ is the reduction of the scintil- 
lation light yield due to escaping electrons. These are 
electrons produced by ionization that thermalize outside 
the Onsager radius and become free from recombination 
even in the absence of an electric field [22J. The effect of 
escaping electrons has been observed for electron recoils 
[22J and has only recently been considered as a possible 
additional factor governing the total scintillation reduc- 
tion for nuclear recoils in LXe [23 J. This is because of the 
surprisingly high ionization yield from nuclear recoils 
HP]. This factor can be expressed in terms of the ratio 
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TABLE I: Xeff and S„ values for the different nuclear recoil energies studied. The third column gives the Xeff values relative to 122 
keV gamma rays. The first error is the statistical uncertainty, and the second error is the systematic uncertainty. 





Er 


-Leff 


Sn 


c 

Sn 


c 

Sn 


c 

Sn 


e 


[keVr] 


at 0.0 kV/cm 


0.73 kV/cm 


l.OkV/cm 


1.5 kV/cm 


4.0 kV/cm 


125 66.7 + 3.3 


1 7S +0.018+0.010 
-0.016-0.009 


0.91 + 0.07 


1.11 + 0.09 


0.88 + 0.06 


_ 


110 57.7 + 3.2 


n 1 09 +0.009+0.004 
-0.009-0.002 


0.95 + 0.05 


0.95 + 0.06 


0.93 + 0.04 


0.93 + 0.06 


95 


46.1 + 4.9 


n 1 CO +0.038+0.010 

W.J-vJO _n noq_n nnq 


0.97 + 0.08 


_ 


0.82 + 0.08 


_ 


80 


34.9 + 2.1 


n 1 QQ +0.042+0.014 


1.33 + 0.26 


_ 


1.30 + 0.25 


_ 


67 


25.7 + 2.0 


n 1 07 +0.025+0.001 
— u.uiy- u.uud 


0.95 + 0.06 


0.91 + 0.12 


0.95 + 0.07 


_ 


58 


19.6 + 2.6 


n 1 +0.046+0.008 
n (T^A nniQ 

— U.UOb— U.UiV 


0.70 + 0.06 


_ 


1.03 + 0.13 


_ 


50 


15.1 + 1.5 


f) 1 70 +0.030+0.019 
w.J.z.o» -0.023-0.014 


0.83 + 0.16 


1.02 + 0.20 


1.01 + 0.15 


1.08 + 0.18 


40 


9.8 + 1.3 


n 1 1 +0.027+0.029 
V.iio _o.o22-0.022 


0.91 + 0.17 


1.64 + 0.50 


0.98 + 0.18 


1.62 + 0.45 


35 


7.6 + 1.2 


n 1 nc +0.028+0.026 
_0.022-0.029 


0.79 + 0.28 


1.06 + 0.30 


0.79 + 0.28 




32 


6.4 + 1.1 


nQ4 +0.027+0.023 
w.wy't _o.o22-0.029 


0.92 + 0.37 


1.25 + 0.45 


0.93 + 0.38 


1.38 + 0.52 


30 


5.7 + 1.1 


D 077 +0.028+0.027 
^•^^^ -0.022-0.026 


1.35 + 0.67 




1.18 + 0.61 




28 


4.9 ± 0.9 


DSH +0.026+0.026 
^•^oo _o.o23-0.032 


1.16 + 0.45 


1.34 + 0.50 


0.87 + 0.35 




25 


3.9 + 0.9 


D 07'^ +0.034+0.018 
u.u/o _o.o25-0.026 


1.19 + 0.52 


1.30 + 0.38 


1.88 + 0.78 





betw^een the initial number of excitons and electron-ion 
pairs a = Nex/M, and the fraction of escape electrons 
over the total electron-ion pairs j8 = Ngsc/M- 



7 esc 



_ Nex + Nj - Nesc _ ^ + 1 - i^NR 



a + 1 - jS 



[122 



(6) 



For this work we select a = 0.06, measured for electron 
recoils in LXe f2T. 24J. is the fraction of escaping 
electrons for 122 keV electron recoils in LXe and is calcu- 
lated to be 0.31 based on the ^^Co data from tlOl and the 
method described in [22 J. ]8nr can be calculated based 
on the nuclear recoil ionization yield reported in 1 10 1 and 
those measured in this work (Figure [14|. The i^nr values 
estimated for the different energies studied in this work 
are given in Figure [l5j 

The last term in the Xeff model is the scintillation light 
quenched by bi-excitonic collisions, qeh as proposed by 
Hitachi 1251 to explain the much lower measured Xeff val- 
ues than are predicted by q^d alone. Bi-excitonic colli- 
sions have the effect of two excitons producing a single 
photon instead of two photons. A recent paper 1261 
extends the study of quenching due to bi-excitonic colli- 
sion by including the varying quenching due to different 
stopping power for different energy recoils, as quantified 
by Birks' Law [27]. Thus, 



1 



dE 
dx 



(7) 



We obtain dE/dx values from SRIM |28| and fit the kg 
parameter to match our Xeff model to the data at 56 keVr, 
finding q^i = 0.65. Hitachi |25 | used qe\ = 0.68 to match 
his model and the data at 60 keVr. 




FIG. 15: Fraction of escape electrons for nuclear recoils, j^NR/ as 
a function of recoil energy used in this work. 



VL DISCUSSION 

The high energy points (56 and 66 keVr) agree with 
previous Xeff measurements. The Xeff measure- 
ments from 20 to 46 keVr suffer from a low differen- 
tial neutron-nuclear elastic scattering cross-section com- 
pared to lower nuclear recoil energies and thus have 
higher backgrounds. Therefore, this results in larger sta- 
tistical errors for measurements in this energy range. 
The detector was designed to minimize multiple scatter- 
ing background for energies below 15 keVr (see Figure 
|9|. For 10 keVr, the Xeff value agrees with previous 
experiments. For energies below 10 keVr, our results 
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are lower than the Chepel f8l results and agree with the 
Aprile [9| values within errors. However, our results 
suggest a decreasing Xeff with decreasing energy and 
not a constant as suggested by [12J. Our results disagree 
with the Xeff curves found by comparin g ca libration 
data and Monte Carlo simulations (FigureJlS) done by 
the XENONIO [12] and ZEPLIN-III [4J collaborations. 
The ionization yield (Figure [l4]) measured in this work is 
in agreement with previous measurements ITOl and the 
XENONIO Monte Carlo analysis \12i . 

The theoretical Xeff models by Lindhard ||2T1 and Hi- 
tachi [25] fail to explain our measurements. The theoret- 
ical model presented in this paper, which includes the 
effect of escape electrons, fits our measurements within 
errors. This Xeff model is a simple one and several cor- 
rections can be done. As pointed out in [29], for LXe 
the Lindhard model is a crude approximation below 15 
keVr, and so the model could be improved below this 
energy. A more general expression for q^sc could also 
be used that includes different values of a for nuclear 
recoils and electron recoils [24] as well as any energy 
dependence of a. 

The Xeff energy dependence measured in this work 
affects the dark matter limits set by LXe detectors. Fig- 
ure [16] shows the XENONIO result [3] (top solid line) 
founa using £^^=0-19, as well as the projected spin- 
independent limits using Xeff=0-19 (solid blue line) and 
the ieff from this work (dotted line). The latter two 
curves are 90% confidence limits based on Feldman- 
Cousins unified approach |30| for a LXe detector with 
a 30,000 kg day exposure with no backgrounds in the 
0.95 to 5.7 keVee energy window. This energy window 
corresponds to 5 to 30 keVr using Xeff=0.19, and to 8.4 
to 39.0 keVr using the measured Xeff in this work. At 
a WIMP mass of 10 GeV/c^, the WIMP nucleon cross- 
section limit is more than an order of magnitude higher 
with the measured Xeff/ compared to the limit found as- 
suming i;eff=0.19. At 50 GeV/c^ the WIMP nucleon cross- 
section limit is only 35% higher and is 20% higher above 
100 GeV/c^. 



collisions and escape electrons that agrees with our re- 
sults. Our Xeff results and model suggest a decreasing 
Xeff with decreasing energy. This result changes the spin 
independent limit as shown in Figure [l6] Although our 
Xeff result significantly changes the cross-section limit at 




WIMP Mass [GeV/cT 



FIG. 16: (Color online) Projected spin independent dark matter 
limits for a LXe detector with 30,000 kg day exposure and 0.95 
to 5.7 keVee energy window. The bottom solid line shows the 
limit with Xeff=0.19 (5 to 30 keVj- window) while the dotted 
line shows the limit calculated with the measured Xeff (8.4 to 
39.0 keVr window). Also plotted are the XENONIO result |3| 
(top solid line) and the regions predicted by |3T1 (red shaded 
region) and i32i (green shaded region). Plot generated using 

m 



low WIMP masses, the limit is only changed by 20% for 
WIMP masses above 100 GeV/c^. 

In addition, we present the results for the nuclear recoil 
ionization yield and calculated the fraction of escape 
electrons for nuclear recoil energies between 4 and 66 
keVr. 



VII. SUMMARY 

This work presents a new £eff measurement for en- 
ergies between 4 and 66 keVr. This measurement was 
done using a single phase (SI signal only) and a dual 
phase detector (SI and S2 signals) to understand the de- 
tector's response at low energies. Each energy measured 
was repeated with at least three different electric fields 
finding no clear dependence of nuclear recoil scintilla- 
tion yield at low energies. We also present a theoreti- 
cal Xeff model including nuclear quenching, bi-excitonic 
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